Ap. J., in press 

A Search for Variability in the Spectral Line Shapes of r Bootis: 
Does this Star Really Have a Planet? 

. . . Artie P. Hatzes and William D. Cochran 

^ I McDonald Observatory, The University of Texas at Austin, Austin, TX 78712 

a^ 

^ ■ ABSTRACT 

O 

CD 

W I An analysis is made of the spectral line shapes of r Bootis using high 

^ I resolution (0.026 A) and high signal-to-noise (S/N ^ 400) data in an effort 

to confirm the planet hypothesis for this star. Changes in the line shape 



are quantified using spectral hne bisectors and line residuals. We detect no 
cn I variations in either of these quantities above the level of the noise in the data. 

One spectral line, Fe I 6213 A, does show a hint of sinusoidal variations in the 
^ I bisector velocity span when phased to the radial velocity period of 3.3 days, but 

this is not seen in the bisectors for two other lines, nor in the line residuals. 
Comparisons of the data to the bisector and residual variations expected for 
nonradial pulsations indicate that we can exclude those sectoral nonradial 
modes having m > 2 and all sectoral modes with k > 1, where k is the ratio of 
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H i the horizontal to vertical velocities for the pulsations. The lack of line shape 

^ I variability and the 469 ms~^ radial velocity amplitude is still consistent with 

nonradial sectoral modes m = 1, and possibly m = 2, but with k ^ 1, which 
k>( ^ is at least 3 orders of magnitude less than the predicted value given the 3.3 

H ' day period of r Bootis. Such low values of k can probably be excluded given 

the lack of photometric variations for this star. Although the measurements 
presented here do not prove, without any doubt, that r Boo has a planetary 
companion, they do add significantly to the increasing body of evidence in favor 
of this hypothesis. 

Subject headings: cool stars - pulsations - planetary companions - variability 



1. Introduction 

The exciting discovery of the planetary companion to 51 Peg (Mayor & Queloz 1995) 
marked a breakthrough in the search for extra-solar planets. After years of unfruitful 
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searches, precise radial velocity (RV) surveys have uncovered several extra-solar planets 
since the 51 Peg discovery (Butler et al. 1997; Cochran et al. 1997; Noyes et al. 1997; 
Butler & Marcy 1996). One of the more unexpected results of these surveys is the fact 
that giant, Jovian-mass planets can exist in very short-period orbits around solar-type 
stars. Although first received with skepticism (primarily because it was so unexpected), 
planet formation theories can now provide a natural explanation for the occurrence of these 
peculiar systems (Boss 1995; Lissauer 1995; Lin et al. 1996). These systems, termed the 
"51-Peg systems" by Butler et al. (1997) are characterized by planets with masses in the 
range of 0.5-4 MjupHer and orbital periods of 3.3-14 days (0.05 - 0.1 AU). Members of this 
class include 51 Peg, r Boo, p^ CNc, and v And. 

Because RV measurements provide only an indirect means of detecting a planetary 
companion, one can never be absolutely certain that the RV signal is due to the reflex 
motion of the star. After all, pulsations, both radial and nonradial, as well as stellar 
surface structure can also produce RV variations which would mimic the signal caused by 
a planetary companion. It is for these reasons that ancillary measurements are needed to 
confirm all planetary discoveries made by RV measurements. Photometric measurements of 
many of the 51 Peg-like stars have established a constant light level to a limit of 0.4 mmag 
(Baliunas et al. 1997). These results provide strong support to the hypothesis that these 
are indeed planetary companions to the 51 Peg-like stars. 

Recent work, however, has questioned whether a planet actually exists around the 
prototypical system, 51 Pegasi. In a controversial paper. Gray (1997) reported changes in 
the spectral line shapes (as measured by the spectral line bisector) of the Fe I 6253 A line 
with the same period as the RV variations in this star. Gray & Hatzes (1997) showed that 
both the changes in the line shapes and the RV variations for 51 Peg could be modeled 
by nonradial pulsations without the need to invoke a planet. Furthermore, the nonradial 
mode would have to be a g-mode oscillation which means that most of the atmospheric 
motion is in the horizontal direction. Since these produce small distortions in the projected 
area of the star, the expected photometric amplitude would be well below detection levels. 
However, the spectral variability for 51 Peg was recently refuted by Hatzes, Cochran, & 
Bakker (1998) using higher quality data; this star most likely has a planet. Even though the 
dispute over the planet of 51 Pegasi has been resolved, it is still worthwhile to confirm the 
other short-period extra-solar planet discoveries and a detailed examination of the spectral 
line shapes may provide the most stringent test for such a confirmation. 

With a semi-amplitude of 469 ms~^ and a period of 3.3 days the "planet" around 
T Boo is the most massive of the 51-Peg systems {m sin i = 3.7 Mjup). Ca II H and 
K measurements imply a rotation period of 3.4 - 4 days for the central star (Baliunas, 
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Sokolof, & Soon 1996; Baliunas et al. 1997), so there is some concern that the RV period 
may be due to rotational modulation. However, the low upper limit (< 1 mmag) to any 
photometric variability of this star argues against surface structure and possibly pulsations 
as a cause for the RV variations (Baliunas et al. 1997). There should also be a concern that 
this star is a nonradial pulsator. Baade & Kjeldsen (1997) found evidence for line profile 
variations in the star e^ Ara that were indicative of an m ~ 4 mode. The spectral type of 
e^ Ara, F6V, is uncomfortably close to that of r Boo (F7V). The RV amplitude for r Boo, 
however, is so large that if it was due to nonradial pulsations then changes in spectral line 
shapes would be large and readily measurable. In this paper we analyze the spectral line 
shapes of r Boo in an effort to confirm or deny the planet hypothesis for this star. 



2. Data Acquisition 

Data were acquired using the 2-D coude spectrograph (Tull et al. 1995) of the McDonald 
Observatory's 2.7-m Harlan J. Smith Telescope. This instrument is a cross-dispersed echelle 
spectrometer used with a Tektronix 2048x2048 CCD detector (24 /im pixels) and it can 
be operated in a "low" resolution {R ~ 60,000) focus with nearly complete wavelength 
coverage in the range 4000 A - 1 /im, or in a "high" resolution focus with more limited 
and non-contiguous wavelength coverage (~ 400 A). The detector was placed at the high 
resolution (Fl) focus for this study. A projected slit width of 2.4 pixels at the detector 
resulted in a working resolving power of 240,000 (resolution = 0.026 A). 

Shortly after commissioning of the 2-d coude we noticed that when using instrument in 
its high resolution mode the theoretical resolving power could not be achieved. Furthermore, 
this limiting resolution seemed to degrade with time. Tests with a pinhole aperture produced 
two images at the camera focus. Most of the time these images fell on a line perpendicular 
to the dispersion which should not affect the instrumental profile (IP) when co-adding all 
CCD rows contained in a spectral order. However, there were instances when these two 
images were not aligned precisely perpendicular to the dispersion, and this would degrade 
the IP. Further tests showed that by masking the top half of the spectrograph collimator we 
could isolate the one pinhole image which remained fixed with time. This also enabled us to 
recover the full resolution of the spectrograph. Work is in progress to isolate and rectify the 
problem (it is suspected that this effect is due to stresses imposed on the cross-dispersing 
prisms in their mount). All observations were made with the top half of the collimator 
masked. This is an investigation of subtle changes in the spectral line shapes which could be 
affected by changes in the IP. Consequently, we chose to sacrifice photons at the expense of 
a stable IP. An examination of the thorium emission lines taken on the nights showed that 
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the rms scatter of the full width to half maximum (FWHM) varied by no more than 10% 
(0.002 A) from observing run to run. Considering that the FWHM of the spectral lines of r 
Boo is about 0.5 A such a small IP variability should not affect the spectral line shapes. 

Table 1 lists the journal of observations which include Julian day of the mid-exposure, 
exposure time, and RV phase. Phases were computed according to the orbital ephemeris of 
Butler et al. (1997): 

Tmax = 2, 450, 235.41 + 3.3128E (1) 

where Tmax represents the phase of maximum radial velocity. 



3. Results 

3.1. Spectral Line Bisectors 

Spectral line bisectors, which consist of the locus of midpoints to horizontal line 
segments extending across the profile, provide sensitive measures of changes in the line 
shapes caused both by surface features (Toner & Gray 1988) or nonradial pulsations (Hatzes 
1996; Gray & Hatzes 1997). However, these types of measurements are less effective for 
stars with rapid rotation since the decreased slope and depths of the lines greatly increase 
the error of the bisector measurement. Also, rapid rotation increases the likelihood for line 
blending. With a rotational velocity of 15 kms^^ r Boo is a less than ideal candidate for 
bisector measurements. Even so, an attempt was made to search for changes in the spectral 
line shapes using line bisectors. 

For a spectral line to be a suitable candidate for bisector measurements it must be 
relatively strong and blend-free. The rapid rotation of r Boo obviously works against these 
criteria. In the spectral region covered by our observations there were only three lines 
suitable for bisector measurements. These are listed in Table 2. 

Prior to computing the individual line bisector the spectra were smoothed with a 
Gaussian having a full width at half maximum of 4 pixels. (Due to the broad lines in r 
Boo we could afford to increase the overall signal-to-noise ratio at the expense of degraded 
resolution.) Individual bisectors (2-4) from a given night were co-added to produce a 
nightly mean. Figures Q and |^ shows the line bisector variations for the Fe I 6213 and Ca 
I 6439 A. (The line bisectors of Fe I 6337 A show comparable variations, or lack thereof, 
and are not shown. The velocity span measurements for this line, however, will be used in 
§3.3). Error bars are shown only for the bisector at phase 0.75 and these represent a typical 
value for the measurements. These were estimated using the noise level as measured in the 
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continuuin and the slope of the flux profiles (Gray 1988). There appear to be no systematic 
variations with phase for the bisectors that is above the noise level. In §3.3 we quantify the 
possible changes in the bisector slopes. 



3.2. Line Residuals 

Because of the line broadening in r Boo, spectral line bisectors may not be the best 
way to search for spectral variability in this star. The rapid rotation, on the other hand, 
may enable us to see directly the distortions in the spectral line profiles due to pulsations. 
It is well known that for stars with rapid rotation both nonradial pulsations and surface 
structure produce distortions that propagate through the spectral line profile. These are 
best seen in stars rotating at vsini > 25 kms~^ , but should still be visible in spectral 
line profiles with the same broadening as r Boo. A convenient means of measuring the 
amplitude of these distortions is through the line residuals (e.g. Kennelly et al. 1992; Baade 
& Kjedsen 1997). 

Line residuals were computed by first co-adding all observations to produce a grand 
mean. Prior to adding the spectra they were all aligned to the rest frame of the star. 
(Again, all observations were smoothed by a 4-pixel Gaussian prior to computing line 
residuals.) Observations on a given night were than combined to produce a nightly mean. 
(The signal-to- noise ratio of these nightly means was ~ 400 per pixel.) The line residuals 
were taken as the difference between the nightly means and the grand mean. Table 3 lists 
the spectral lines that were used to search for residual variations. Figure ^ and Figure ^ 
are representative residuals for the four spectral lines in the wavelength interval 6108 A - 
6117A and the Fe I line 6224 A. There appears to be no phase-dependent distortions across 
the spectral line profile and any distortions which may be present are at the level of the 
noise (as indicated by the scatter of the residuals in the continuum). 



3.3. Limits on Nonradial Pulsations 

At the present time, nonradial pulsations seems to be the only hypothesis, other than 
a planetary companion, that may provide an explanation for the RV variability of r Boo. 
Even though the estimated rotational period for r Boo is near that of the RV period, surface 
features (spots, etc) can be excluded as a possible hypothesis not only on the basis of the 
lack of photometric variations, but also from the lack of spectral distortions as evidenced 
in our bisector and line residual measurements. It is difficult for any surface feature to 
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account for the large RV amplitude in r Boo without producing measurable changes in the 
line shapes. 

It is possible, however, for nonradial pulsations to produce significant RV variability 
without large variations in the spectral line shapes (Hatzes 1996). The question is, 
whether these can produce the large RV amplitude seen in r Boo without measurable 
line distortions. In this section we place limits on any possible nonradial pulsation mode 
which may be present on r Boo. We considered all sectoral nonradial modes that were 
capable of producing the observed RV amplitude. Furthermore, the shapes of the RV curves 
from nonradial pulsations are nearly sinusoidal, consistent with the nearly zero orbital 
eccentricity. 

The expected line bisector and line residual variations from nonradial pulsations were 
calculated using the prescription outlined in Hatzes (1996). Nonradial pulsations can be 
described in terms of spherical harmonics Y^^m where d, is the degree and m the azimuthal 
order. Only sectoral modes (£ = m) were considered as these should produce the largest 
integrated RV amplitude. The pulsational amplitude for each mode was chosen such that 
the integrated RV amplitude from the oscillations was 469 ms~^ , i.e. the observed RV 
amplitude in r Boo. 

In calculating the local velocity field on the stellar surface due to nonradial pulsations 
one needs to know the ratio of the horizontal to vertical velocity amplitudes of the 
atmospheric motions. Normally this is not a free parameter, but can be calculated by 

k = (^^ (2) 

V0.116/ ^ ' 

where Q is the pulsation constant = Ii{p/ PqY^"^ . (11 is the pulsation period in days, p is the 
mean density of the star, and p© is the mean density of the Sun). For r Boo the appropriate 
value of k is ~ 1000. However, there is a well know "k-parameter problem" where in 
modeling the line profile variations of known nonradial pulsators, the profiles can be better 
fit with low values of k (~ 0) even though the pulsation period implies the need for much 
larger values of k (Lee & Saio 1990). Consequently, k was treated as a free parameter. 

A convenient means of quantifying a line bisector is through its velocity span. Due to 
the large errors in the bisector data we chose to measure this velocity span using a linear 
regression of all bisector points within two arbitrarily chosen flux points. As can be seen 
from Figs. 1 and 2 measuring the span using only two bisector points can result in grossly 
different velocity span measurements depending on which end points are chosen. The slope 
of the linear regression through all bisector points is less sensitive to the choice of end 
points. Comparisons of this span to those predicted by nonradial pulsations were done in a 
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consistent manner, i.e. the bisector from a synthetic spectral hne profile having nonradial 
pulsations was computed and the velocity span was calculated using a linear regression 
through all bisector points within the specified flux values used in the data bisectors. The 
synthetic data, like the real data, were also smoothed by a Gaussian having a FWHM of 4 
pixels. 

Figures |^ - show the bisector velocity span measurements for the Fe I 6213 
A, Fe I 6337 A, and the Ca I 6439 A lines, respectively. The averaged velocity span 
measurements are tabulated in Table |[ The error was estimated by taking the standard 
deviation, cr, of the individual span measurements clustered around phase 0.45. This ±a 
was plotted as the error on each of the data points. Superimposed on these figures are the 
predicted velocity span variations from nonradial modes having k=10, m=l,2 and k=1.2 (a 
canonical value for g-mode oscillations), m=l,2. Use of the predicted k value {k ^ 1000) 
for r Boo resulted in span variations much larger than those shown in the figures. This was 
also true for sectoral modes with m > 2, regardless of the fc-value. We can thus exclude 
those pulsation modes with m > 2 and all modes with k ^ 10. Only the m = 1 with k = 
1.2 would produce bisector span variations which would be hidden by the noise in our data. 
The presence oi an m = 2, k = 1.2 mode is marginally possible and, in fact, the Fe I 6213 A 
bisector spans do show a hint of sinusoidal variations with an amplitude comparable to that 
expected for a.n m = 2 (albeit out-of-phase). However, we deem this to be merely noise as 
there are no obvious sinusoidal variations in the span measurements for the other two lines. 

Even though the phased bisector velocity span measurements show, for the most 
part, no obvious sinusoidal variations, a periodogram analysis was still performed these 
measurements (Table ^. There was no statistically significant power at the RV period 
of 3.31 days in any of the periodograms. The false alarm probability for any power at 
the RV period was 0.4 for the Ca I 6439 A line and 0.2 for the Fe I 6337 A line. Power 
did appear near the appropriate frequency in the periodogram for the Fe I 6213 A span 
measurements with a moderately low false alarm probability of 0.02. Again, we do not 
deem this significant based on the lack of power at the appropriate period in the other 
periodograms. (We should point out that the best line for bisector measurements, the Ca I 
6439 A, showed the highest false alarm probability.) 

Nonradial pulsations should also be visible in the spectral line residuals. Figure |^ 
shows the percent change in the model spectral line residuals for an m=2, k=10 mode. The 
distortions have an amplitude of a few percent, much larger than the observed distortions 
in Figure ^ and ^ (The amplitude of the distortions for the m = 1, k = 1.2 mode are about 
0.5%.) Figure ^ shows the standard deviation of the predicted spectral line residuals for 
various pulsation modes as a function of RV phase. The standard deviation was computed 



over the full width of the spectral line. The crosses represent the standard deviation of the 
line residuals of the actual data, again computed over the full width of the spectral line. 
Individual measurements are shown as crosses while the solid dots represent a mean value 
for all seven spectral lines that were examined. Part (if not all) of the rms scatter that was 
computed for the real data may be due to photon noise. In order to correct for this, the a 
over an equivalent wavelength region in the continuum was calculated and subtracted in 
quadrature from the a computed across the spectral line. 

It is not strictly correct to subtract the a in the continuum since the signal-to-noise 
ratio in the spectral line decreases from the continuum value. For line depths of 10-15% 
and SjN ^ 150, this effect should result in a value a ~ 0.0015 for the residuals across 
the line profile. This is comparable to the mean level of a for the line residuals as seen in 
Figure |^. We therefore conclude that the non-zero mean a for the line residuals can be 
explained entirely by noise and is not due to intrinsic variability of the star. 

In comparing the variation of a for the line residuals with the pulsation models we can 
once again exclude all pulsations modes with m > 2 and those with k ^ 10. Of the modes 
having the smaller horizontal scaler, only an ?7i = 1 is consistent with the apparent lack of 
line shape variations. An m = 2 can be hidden in the residual variations for an individual 
line, but it seems to be excluded when using the residual variations averaged over all lines. 



4. Discussion 

Our analysis of the spectral line shapes of r Boo, both through line bisectors and line 
residuals, indicate no variability above the noise of the measurements. This lack of spectral 
variability is consistent with recent findings for this star by Brown et al. (1998). If r Boo 
is really a pulsating star, then the only pulsation mode it could possibly have is an ?7i = 1, 
but with a low horizontal scaler {k ^ 1), a value several orders of magnitude less than the 
predicted value. This mode and possibly an m = 2 (although less likely) could account for 
the observed RV variations in this star and yet show no measurable changes in the spectral 
line shapes. 

One would think that these modes could be excluded on the basis of the low k values, 
but this may not be the case in light of the fc-parameter problem discussed earlier. A 
calculation of pulsations in a real stellar atmosphere along with the temperature variations 
may result in a much lower "observed" A;- value. 

If the low k values implied by this study do indeed represent a true ratio of the 
horizontal to vertical velocities, then the pulsational amplitude can be used to estimate the 
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photometric amplitude for the oscillations. For both the m = 1 and m = 2 modes {k = 1.2) 
the pulsational amplitude is about 400 ms^^ . Buta & Smith (1979) published expressions 
for estimating the photometric amplitude due to geometrical effects. These result in AV ~ 
0.03 mag for the m = 2 mode and AV ~ mag for the m = 1 mode. 

There may, of course, also be photometric variations due to temperature effects and 
for very long periods these may dominate the light curve (Buta & Smith 1979). The 
temperature variations can be estimated using 



1 ^ UO Ve 



where \/ad is the adiabatic temperature gradient, Vt is the rotational frequency, u is the 
pulsational frequency, Vp is the pulsational amplitude and Ve is the rotational velocity (Lee 
& Saio 1990). For r Boo, uo k, Vt^ Vp = 0.4 kms""*^ , and v^ = 15 kms~^ . For the m {= i) = 
1 mode (5T ?^ 170 K, and for the m = 2 mode 6T ^ 300 K. These variations result in AV ~ 
0.1 and 0.2 mag, respectively, for the two modes, assuming that these temperature changes 
represent disk integrated quantities. This is a reasonable assumption given the low degree 
of the modes. The lack of photometric variations seems to exclude such large temperature 
variations. However, as pointed out by Buta & Smith, the temperature variations for very 
long period modes (i.e. much greater than the period of the fundamental radial mode) 
are difficult to calculate and these may have to include both non-adiabatic and non-linear 
effects. Consequently, these approximate estimates, derived from linear theory, may not be 
appropriate for the hypothetical pulsation mode in r Boo. Furthermore, Buta & Smith 
state that the light curve from temperature variations may be in-phase or out-of-phase with 
that due to the geometric variations. One could imagine a scheme where the geometric 
variations nearly cancel those due to temperature variations producing a net zero change 
in the light level. These may also conspire to produce small changes in the spectral line 
shapes, even if k is large. Clearly, a detailed analysis of pulsations in a stellar atmosphere 
is required to obtain better estimates of any predicted light curve. At face value, however, 
the photometry does seem to exclude the possibility of such low degree nonradial modes. 

Of course there is no reason to expect that the hypothetical pulsation modes are 
spheroidal in nature. One possibility is toroidal, or r-modes (Saio 1982). Because the 
atmospheric motions for these modes are purely horizontal, they should not be accompanied 
by photometric variations. It is not clear, however, if these modes could produce appreciable 
RV variations. For a star with intermediate inclination r-modes can produce spectral line 
distortions (Unno et al. 1989), but these are not seen in our spectral data. 

Astrophysics, more than most sciences, has its level of uncertainty. This is primarily 
because the objects of our study are distant and cannot be scrutinized in laboratories. 
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or up-close in exquisite detail. The best we can do is collect high quality data and draw 
the "most probable" hypothesis that is consistent with that data. There will always be 
some probability, no matter how small, that the conclusions that are drawn are incorrect 
due to our ignorance of certain aspects of natural phenomena. This is especially true 
for extra-solar planet discoveries made with RV measurements since these are indirect 
detections and some unknown process may be the cause of any observed variability. One 
thing is certain, if we had found spectral variability in the line shapes of r Boo with the 
RV period, then this would almost certainly have killed the planet hypothesis. Although 
our analysis shows no spectral variability in r Boo it does not prove beyond any doubt 
that this star has a planetary companion. However, our measurements certainly add to the 
considerable evidence weighing in favor of the planet hypothesis. It is still the simplest and 
most likely explanation for the RV variability in this star. 

Of the 5 short-period extra-solar planets, two have been confirmed using analyses of 
the spectral line shapes. Although it seems highly likely that all the short-period systems 
have planets, it is still worthwhile to search for any spectral variability in these other stars. 
We are in the process of obtaining high resolution measurements for the other short-period 
planets in order to confirm the planet hypotheses for these stars. 

This work was supported by NASA grant NAGW3990 and NAG5-4384. 
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Table 1. Journal of Observations 



JD 


Phase 


S/N 


JD 


Phase 


S/N 


2450481.855 


0.392 


110 


2450530.900 


0.196 


135 


2450481.890 


0.402 


125 


2450530.938 


0.208 


120 


2450481.898 


0.405 


135 


2450530.973 


0.218 


130 


2450481.933 


0.415 


130 


2450531.722 


0.445 


175 


2450482.939 


0.719 


125 


2450531.730 


0.447 


160 


2450482.977 


0.730 


130 


2450531.768 


0.458 


140 


2450482.982 


0.732 


120 


2450531.803 


0.469 


150 


2450528.816 


0.567 


210 


2450531.939 


0.510 


120 


2450528.850 


0.578 


170 


2450531.974 


0.521 


130 


2450528.858 


0.580 


180 


2450584.732 


0.446 


170 


2450529.767 


0.858 


150 


2450584.761 


0.455 


185 


2450529.773 


0.856 


160 


2450585.704 


0.740 


175 


2450529.809 


0.867 


170 


2450585.766 


0.758 


200 


2450529.822 


0.871 


190 


2450585.774 


0.761 


160 


2450530.844 


0.180 


140 


2450585.808 


0.771 


150 


2450530.853 


0.182 


170 


2450586.728 


0.049 


230 


2450530.892 


0.194 


150 


2450586.763 


0.059 


170 



Table 2. Spectral Lines for Bisector 
Measurements 



Wavelength 

A 


Species 


EW 
(mA) 


6213.3 


Fel 


75 


6336.9 


Fel 


90 


6439.1 


Cal 


160 



Table 3. Spectral Lines for Residual 

Measurements 



Wavelength 

A 


Species 


EW 
(mA) 


6108.1 


Nil 


55 


6110.5 


? 


30 


6113.3 


Fell 


40 


6116.2 


Nil 


70 


6224.0 


Nil 


25 


6335.4 


Fel 


95 


6336.9 


Pel 


55 



Table 4. Velocity Span Measurements 



JD 


Phase 


'S'6213 


'S'ess? 


'S'6439 
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Fig. 1. — The spectral line bisectors as a function of phase for Fe I 6213. The error bars 
shown at phase 0.75 represent typical errors. 
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Fig. 2. — The spectral line bisectors as a function of phase for Ca 1 6439. 
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Fig. 3. — (Top) The spectral line residuals as a function of phase for the 6107 - 6117 A 
spectral region of r Boo. (Bottom) The grand mean spectrum in the same wavelength 
interval. 
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Fig. 4. — The spectral line residuals for Fe I 6224 A (top) and the grand mean spectrum 
(bottom). 
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Fig. 5. — Bisector velocity span for Fe I 6213 A. The lines show the predicted variations for 
various pulsation modes. 
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Fig. 6. — Bisector velocity span for Fe I 6337 A compared to the predicted variations from 
nonradial pulsations. 
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Fig. 7. — Bisector velocity span for Ca I 6439 A compared to the predicted variations from 
nonradial pulsations. 
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Fig. 8. — Predicted spectral line residuals as function of phase for an m=2 nonradial mode 
with k = 10. 
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Fig. 9. — The crosses represent the standard deviation, a, of the hne residuals across the 
spectral line profile for the seven spectral features that were examined. The solid points 
represent the average of all seven lines. The lines represent the predicted o" as a function of 
phase for various nonradial modes. 



